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We study the consequences of negative differential electron mobility in insulated gate field effect 
transistors (FEFS) using the field model. We show that, in contrast to the case of the monotonic 
velocity saturation model, the field distributions in a short-channel FET may be described by 
the gradual channel approximation even for high drain-to-source voltages. The current-voltage 
dependence of the short-channel FET should have a branch with a negative slope. The FET exhibits 
a negative differential resistance and may show convective or absolute instability, depending on the 
applied voltages. The fluctuation growth is governed by the diffusion law with a negative effective 
diffusion coefficient. 
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INTRODUCTION 

Since the electric fields in modern short-channel field 
effect transistors (FETs) may be very high, effects of elec- 
tron velocity saturation are of importance. Studies im- 
plying monotonic saturation of the electron velocity with 
increasing electric field [l]-[4] have shown that in short- 
channel FETS the drain saturation current j sat and the 
drain-to-source voltage U sa t at which the current satu- 
ration takes place are substantially smaller than the sat- 
uration current and voltage predicted by the constant 
electron mobility model. When the drain-to-source volt- 
age Ud is smaller than U sa t , the drain current increases 
with increasing Ud and may be described by the grad- 
ual channel approximation [4] . For Ud > U sa t (in the 
current saturation region), the gradual channel approxi- 
mation fails because of too sharp variation of the electric 
field along the channel. At this takes place, the drain 
current is supposed to be constant. 

Recently, Dmitriev, Kachorovskii, Shur and Stroscio 
[5] have shown that two-dimensional electron gas in sil- 
icon or germanium should exhibit a negative differen- 
tial mobility caused by the electron runaway. Appar- 
ently, the above consideration implying monotonic elec- 
tron velocity saturation is not sufficient to describe FETs, 
which may exhibit non-monotonic velocity-field depen- 
dence with negative differential electron mobility. This 
may be the case in Si or Gc MOSFETs or FET based on 
compound semiconductors, such as GaAs and InP. 

In metal semiconductor FETs (MESFETs), the conse- 
quences of the negative differential mobility arc similar 
to ones in bulk devices resulting in the Gunn effect [6] . In 
particular, the stationary Gunn domain which may be de- 
scribed in a three-dimensional model plays an important 
role in MESFETs [7]. In contrast, the consequences of 
the negative differential mobility should differ essentially 
in two-dimensional case such as metal-oxide semiconduc- 
tor FETS (MOS-FETS), metal-insulator semiconductor 
FETS (MISFETS) and heterostructure FETS (HFETS). 



Indeed, the distribution of electric field caused by a 
charged fluctuation is quite different in 2D layer and 
also depends on the induced charge in the conductive 
gate electrode. Due to this fact, the relaxation of a 
charged fluctuation is governed by the diffusion law [8], 
in contrast to conventional exponential law for three- 
dimensional Maxwellian relaxation. However, no ana- 
lytical theory describing the consequences of the negative 
differential mobility in two-dimensional FET channel has 
yet been developed. 

In this paper, we study the consequences of the neg- 
ative differential electron mobility in the insulated gate 
FETS, such as MOSFETS, MISFETS and HFETS, using 
the field model and the gradual channel approximation. 
In these devices, the channel is isolated from the gate by 
an insulator( a wide band gap semiconductor in HFET 
). We shall not consider here MESFETS and junction 
FETS (JFETS) that have differing design and require 
other models. 

We have found that the FET may exhibit convective 
or absolute instability, depending on the applied voltages. 
The growth of small electric fluctuations is governed by 
the diffusion law with a negative effective diffusion coef- 
ficient. We have also shown that, in contrast to results 
following from the monotonic velocity saturation model, 
the current-voltage dependence of the short-channel FET 
should have a branch described by the gradual chan- 
nel approximation even at high drain-to-source voltages. 
This branch has a negative slope, i.e. the FET exhibits 
a negative differential resistance. 



BASIC EQUATIONS 

The surface electron concentration n in the FET chan- 
nel is determined by the conventional charge control 
model equation [4] 
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where C is the surface gate capacitance, Uac( x ) is the 
gate-to-channel voltage, Ut is the threshold voltage. 
Eq.(l) is valid in the gradual channel approximation 
when the characteristic scale of variations of U and n 
is much greater than the gate-to-channel separation d. 

In view of Eq.(l), the continuity equation may be writ- 
ten as 
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Here, V(E) is the electron drift velocity which depends 
on the longitudinal electric field E in accordance with 
the local field model. The V{E) dependence in the 
FET channel may differ quantitatively from that in the 
bulk semiconductor but have a similar form with a neg- 
ative slope for E greater than the threshold field E t (see 
Fig.l,a). The dynamics of E distribution is thus de- 
scribed by Eq.(2) with the use of Fig.l,a. 



STEADY STATES 

In a steady state, Eq.(2) becomes U = j/CV where j 
is the drain current per unit gate width. Differentiating 
this expression, we obtain 
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where fx d is the differential electron mobility. The bound- 
ary values of E, U and V at the source (x = 0): E ,Uo 
and Vq are determined by the condition 



V(Eo) = V = 
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Uo = U G - U T , 
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where Uq = Ugc{0) is the gate voltage, which is assumed 
to be fixed. Eqs.(3) and (4) allow one to calculate the 
E distributions for certain j if the V{E) dependence is 
given in the form of an analytical approximation. One 
can however analyze their solutions qualitatively without 
calculations by simply looking at the (dE/dx, E) diagram 
( see Fig.2,b) which results from Eq.(3) and Fig.l,a. It 
will be seen that E increases along the channel when 
E < E t and drops when E > E t . Note that, at E > E t , 
dE/dx increases with E faster than E. Therefore, E 
drops from infinity at a finite distance x m . The general 
form of a solution of Eq.(3) is presented in Fig. 2. 

As it follows from Eqs.(3),(4), the characteristic scale 
of field variation a is given by a = (Vo /V t )(Uo / E t ), where 
Vt = V(E t ) (in particular, x m is of the order of a). As- 
suming Uq ~ IV, E t ~ 3 kV/cm, Vo ~ V t , we have 
a ~ 3mkm. Note that the gradual channel approxima- 
tion used here is valid when o>d. 

The actual field profile in the channel is represented 
by the section of the E(x) curve( see Fig. 2 ) which ex- 
tends to the right from the point where E = Eq for a 
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FIG. 1: Schematic V(E) dependence (a), (dE/dx, E) dia- 
gram (b) and current-voltage characteristics ](Ud) for differ- 
ent channel lengths; broken lines correspond to the possible 
current saturation regime which is not described by the grad- 
ual channel approximation (c). 



distance of the channel length L. As it can be seen from 
Fig.l,a, if V s < Vq < V t , two values of E satisfy Eq.(4): 
Eqi < E t and E 2 > E t . Then, two steady states de- 
termined by boundary fields Eoi and E02 may exist in 
a FET with given drain current j: a positive state with 
E < Et, dE/dx > 0,fid > and a negative state with 
E > E t , dE/dx < 0, /id < (see inset in FigEJ. These 
two states differ in the value of the drain-to-sourcc volt- 

L 

age Ud = J Edx. It follows from Fig.l,a and Fig. 2 that 


the negative state can exist if the channel length L is 
short enough: L < L m . The maximum length L m de- 
pends on j: L m — > when Vq = j/(CUo) — > Vt, and 
L m — > x m — a when V — ► V s . 

The current voltage characteristics j(Uo) following 
from this analysis are given in Fig.l,c. For short enough 
channel length L < x m , the j(Ud) dependence has two 
branches: I and II, corresponding to the above mentioned 
positive and negative states, respectively (curves I and 
2). The type II branches existing for high drain-to-source 
voltages Ud > E t L show a negative differential resis- 
tance. In the voltage interval near Ud ~ E t L, no steady 
states exist in the considered model. Actually, the cur- 
rent saturation regime similar to the one considered in 
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above, linearizing Eq.(2) and using Eq.(3), we obtain 
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FIG. 2: E(x) dependence representing the general solution 
of Eq.(3). Inset: steady field profiles along the channel with 
given drain current. 



the monotonic velocity saturation model [3] should take 
place, which is not de-scribed by the gradual channel ap- 
proximation. It is possible that states corresponding to 
the current saturation regime may also exist at higher 
Ud in the region where branch II lies (see broken lines in 
Fig.l,c). Then the FET may presumably show branch II 
or the current saturation branch, depending on the cir- 
cuit condition, applied voltages and the previous state 
of the FET. The question as to which of these branches 
should be observed requires three-dimensional computer 
simulations and is beyond the scope of this paper. 

For very short channel lengths d <C L <C a, the j(Ujj) 
characteristic is close to the one determined by the V(E) 
dependence since the field in the channel is close to its 
boundary value Eqi or E 02 (curve 1, j = CUqv(Ud/L)). 
For L > x m , branch II is absent and current saturation 
should occur at high Ud (curve 3) as is the case in the 
conventional model. 



STABILITY OF THE STEADY STATES 

Let us consider the temporal behavior of small electric 
fluctuations in the channel. Adding small corrections E 1 
and U' to the steady distributions of E and U described 
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D = U^ d , S = V + fi d E-U 
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The parameters r, D, S and \Xd vary with the charac- 
teristic scale a. 

It can be shown that in the Fourier spectrum of fluctu- 
ations, only components with high k arc of importance: 
a -1 <C k <C d^ 1 . Then, they behave as exp(ikx — iujt) 
where ui varies slowly over the length fc -1 . Putting 
this in Eq.(5) and noting that rDk 2 3> 1, we find: 
oj = kS—iDk 2 . Note that this expression and Eq.(5) look 
like equations describing diffusion and drift processes, D 
and S being effective diffusion coefficient and drift veloc- 
ity, slowly varying with coordinate x. Thus, fluctuations 
behave in accordance with the diffusion law. It can be 
seen that the negative states (branch II) are unstable 
since the effective diffusion coefficient D = U \x& is neg- 
ative. Then fluctuations grow with an increment |_D|/c 2 . 
The fast growing components with k ~ d" 1 arc most 
essential because the gradual channel approximation is 
invalid for bigger k. The positive states are stable and 
fluctuations are damped out there since D > 0. As we 
lave shown in [7] , a similar diffusion relaxation with con- 
stant D = fj,Uo takes place in the FET at Ud = 0, fi being 
the low field mobility. 

The consequences of this convective instability depend 
on whether fluctuations have a chance to grow essen- 
tially during their non-uniform motion along the channel. 

L 

Their growth is governed by exp(g),g ~ k 2 J \D/S\dx, 

o 

k ~ d . Under the condition g < 1, analogous to the 
well-known Kroemer criterion for the Gunn diode, the 
FET should presumably act as a stable amplifier, and at 
g > 1 as a generator. Note that, as it can be seen from 
Eqs.(3) and (5), the velocity S becomes zero at a certain 
value E > E t . If the field in the channel reaches this 
value, absolute instability occurs i.e. fluctuations grow 
without motion near the corresponding point x in the 
channel. Detailed study of the instability criterion, re- 
quiring a numerical calculation, is beyond the scope of 
this paper. 
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